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ABSTRACT 

Observations of the Seyfert 2 and starburst galaxy NGC 5135 with the Chandra X-ray Observatory 
demonstrate that both of these phenomena contribute significantly to its X-ray emission. We spatially 
isolate the active galactic nucleus (AGN) and demonstrate that it is entirely obscured by column 
density Nh > lO^^'cm"^, detectable in the Chandra bandpass only as a strongly reprocessed, weak 
continuum and a prominent iron Ka emission line with equivalent width of 2.4 keV. Most of the soft X- 
ray emission, both near the AGN and extending over several-kpc spatial scales, is collisionally-excited 
plasma. We attribute this thermal emission to stellar processes. The AGN dominates the X-ray 
emission only at energies above 4 keV. In the spectral energy distribution that extends to far-infrared 
wavelengths, nearly all of the emergent luminosity below 10 keV is attributable to star formation, not 
the AGN. 

Subject headings: Galaxies: individual (NGC 5135) — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

Active galactic nuclei and star formation are fun- 
damentally related. In their quiescent states, black 
hole and stellar components are observed to be re- 
lated in individual galaxies in the rough proportional- 
ity between the black hole mass and the velocity dis- 

I persion of the stellar sp heroid ijGebhardt et al.l [20001: 

, iFerrarese fc Merriti]l200rl) . suggesting a connection over 
cosmological timescales. Also, Kauffmann et al.l 1)20031) 
find a strong connection between active galactic nucleus 
(AGN) luminosity and the age of the stellar popula- 
tion of the host galaxy for AGNs in the Sloan Digital 

[ Sky Survey. Locally, AGNs and strong star formation 
are correlated, with roughly half of well-selected sam- 
ples of Seyfert 2 galaxies exhibiting compact circum- 
nuclear starbursts llCid Fernandez ; et aL ll2f)0Tl and ref- 
erences therein). These Seyfert 2/star burst composite 
galaxies contain genuine AGNs, which accretion onto the 

' central black hole powers, but the stellar processes can 
also become energetically significant. 

The local composite galaxies afford detailed exam- 
ination for discrimination of their AGN and stellar 
contributions and accurate measurement of these com- 
bined effects. The detailed studies serve as an impor- 
tant foundation for understanding ultraluminous infrared 
galaxies, whose bolo metric luminosities exceed lO-'^^L© 
iSanders fc MirabelllTool . Among this class, the prob- 
lem is to identify either an AGN or star formation as 
the primary energy source. If present, an AGN may be 
deeply buried, obscured both on small scales associated 
directly with the central engine an d on the large scales 
of star formation ijPtak et al.ll2003() . While the compos- 
ite galaxies are not the direct analogues of these com- 
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plex higher- luminosity cases, they serve as basic build- 
ing blocks, demonstrating the primary consequences of 
combining starbursts with AGNs in highly obscured cir- 
cumstances. 

Distinguishing these X-ray emission processes is also 
important in the study of the cosmic X-ray background 
(XRB). AGN can fun damentally produce the XRB (e.g., 
iGiacconi et al. 11200 ll|l . but they require a distribution of 
redshi ft and absorption to m atch its observed spectrum 
(e.g., iSetti fc Wolti er"1989^. The outstanding difliculty 
is to identify the highly-obscured populations, which are 
not always evident at energies below 10 keV. The Seyfert 
2/starburst composite galaxies cannot directly account 
for the XRB, but they tend to be s trongly obscured 
ijLevenson. Weaver, fc Heckmanll2001a|) . and they illus- 
trate the systematic effects of observing AGNs in the 
presence of starbursts, even at X-ray energies. 

With high spatial resolution and simultaneous spec- 
troscopic data from the Chandra X-ray Observatory, 
we analyze here the X-ray emission due to both AGN 
and starburst components of the Seyfert 2/starburst 
composite galaxy NGC 5135. This galaxy is opti- 
cally classified as a Seyfert 2 on the basis of e mis- 
sion line ratios l)PhilliDS. Charles, fc Baldwin! 1198^ . It 
also contains a powerful starburst within 200 pc of 
the nucleus that is most evident at ultraviolet ener- 
gies l|Gonzalez Delgado et al.llT998l) . NGC 5135 is rel- 
atively nearby, at distance of 59 Mpc (assuming Hq = 
70 kms^^ Mpc~^), so 1" corresponds to 285 pc. In X- 
rays, we expect to find the characteristic unresolved 
non-thermal con tinuum from accre tion onto the cen- 
tral black hole ijTurner et alJ Il997j) . We also distin- 
guish the dominant X-ray signature of a starburst: spa- 
tially extended thermal emission, due to individual su- 
pernovae and stellar winds, which collectively may pro- 
duce a "superwind" outfiow that escapes the galaxy 
ijDahlem. Weaver, fc Hec kman 19981. With the simulta- 
neous spatial and spectral distinction of these emission 
processes, we can more accurately measure them inde- 
pendently. Relating the X-ray data to multi-wavelength 
observations of this galaxy as a whole, we draw general 
conclusions about identifying multiple energy sources 
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in low-resolution observations of galaxies that contain 
AGNs. 

2. OBSERVATIONS AND DATA REDUCTION 

The Chandra Advanced CCD Imaging Spectrometer 
(ACIS) observed NGC 5135 on the back-illuminated S3 
detector on 2000 September 4. (See Wcisskopf et alJ200nl 
for more information on Chandra.) We reprocessed all 
data from original Level 1 event files removing the spec- 
tral and spatial randomization that is included in stan- 
dard processing with Chandra Interactive Analysis of Ob- 
servations (CIAO) software. We applied the current gain 
corrections for the S3 CCD (from 2001 July 31), and 
we included only good events that do not lie on node 
boundaries, where discrimination of cosmic rays is diffi- 
cult. We examined the lightcurves of background regions 
and found no significant flares, so we used all data within 
the standard good-time intervals, for a total exposure of 
29.3 ks. 

3. IMAGE ANALYSIS 

These Chandra data reveal the complexity of NGC 
5135's X-ray emission, as we illustrate with the total 
0.4-8.0 keV image (Figure^). The raw (O'.'S) pixels have 
been smoothed by a Gaussian of FWHM=1'.'5, and both 
the image intensity and contours are scaled logarithmi- 
cally. The active nucleus is not the only source of X-rays. 
Rather, the image reveals diffuse emission, off-nuclear 
point sources, and, most surprisingly, two bright cen- 
tral sources, separated by 2". As we demonstrate below, 
only the northern of these two is the active nucleus. The 
highly-smoothed version of the image (Figure shows 
that the soft emission extends over a scale of approxi- 
mately 50" X 30" (14 kpc X 9 kpc), tracing the large-scale 
bar in the disk of NGC 5135. 
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Fig. 1. — Chandra X-ray image of NGC 5135. The broad- 
band (0.4—8 keV) image has been smoothed by a Gaussian of 
FWHM=1'.'5. The image is scaled logarithmically. The overlaid 
intensity contours begin 60" above the background level and in- 
crease by factors of two. The AGN is the northern central source, 
while the starburst is responsible for the southern nuclear concen- 
tration and diffuse emission. Several additional sources are also 
apparent. 



We consider three broad energy regimes separately: 
0.4-1.0, 1.0-4.0, and 4.0-8.0 keV, which we refer to as 
the soft, medium, and hard images. In soft X-rays, (Fig- 
ure the galaxy appears very similar to the total im- 
age. The very extended emission is soft, and both central 
sources are evident. Some of the additional sources to- 
ward the east are less prominent or do not appear at 
lower energies. Some of the extended emission is still de- 
tected in the medium X-ray (1-4 keV; Figure 2J image, 
and the additional sources are more prominent. In hard 
X-rays (4-8 keV; Figure [SI, the northern central source 
is strongest. The southern source contributes weakly, 
but is still evident, and two of the eastern sources are 
significant. 
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Fig. 2. — Digitized Sky Survey red image with 0.4-8.0 keV 
X-ray contours overlaid. The X-ray data have been smoothed by 
a Gaussian of FWHM=5" to illustrate the correlation of X-ray 
and optical emission on large scales. The optical image is scaled 
linearly, and the X-ray contours are scaled logarithmically. 



We compare the two central sources with Chandra^s 
point-spread function (PSF) in the three separate en- 
ergy bands. In each of these three bands, we model the 
emission profile as a constant plus two two-dimensional 
Gaussians having circular cross-sections. In all cases, 
the Gaussians are significantly broader than a similar 
two-dimensional fit to the PSF at the appropriate en- 
ergy and detector position modeled with the CIAO task 
mkpsf. We conclude that the two central sources are 
resolved in all three broad bands. Therefore, neither of 
the broad-band emission peaks is due to the unresolvable 
AGN alone. 

We do, however, morphologically identify the AGN in 
the Fe Ka emission of NGC 5135. This fluorescent line 
radiation arises near the strong AGN continuum and 
advantageously avoids confusion with the soft, thermal 
starburst. In this case, the Fe line is prominent above the 
continuum that either the AGN or other sources could 
produce. Isolating the emission in the 5.5-7.5 keV band- 
pass and subtracting the continuum measured from 4 to 
5.5 keV produces an image in Fe Ka. This image con- 
tains an unresolved point source, which is located 0"5 
north of the northern emission peak that is measured in 
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Fig. 3.— Chandra soft X-ray (0.4-1 keV) image of NGC 5135, 
smoothed by a Gaussian of FWHM=1'.'5. The image is scaled 
logarithmically. The overlaid intensity contours begin 8cr above 
the background level and increase by factors of two. 



Fig. 4. — Chandra medium X-ray (1-4 keV) image of NGC 
5135, smoothed by a Gaussian of FWHM=1'.'5. The image is scaled 
logarithmically. The overlaid intensity contours begin 6(t above the 
background level and increase by factors of two. 



the total band image. Thus, the Fe Ka image specifically 
locates the AGN on the northern side of the broad peak 
that appears in the total band image, with both the AGN 
and other sources contributing to the spatially-extended 
broad-band peak. At the distance of NGC 5135, the 
point-like morphology constrains the physical scale of the 
Fe Ka emission to i? < 285 pc. 

4. SPECTROSCOPY 

We extracted spectra from several interesting regions. 
In all cases, we also measured the local background in 
nearby source-free regions and subtracted it. Over the 
course of the Chandra mission, the soft X-ray sensitivity 
has diminished, likely the result of build-up of material 
on the detector. We have used the "ACISABS" model of 
G. Chatras and K. Getman^ to create ancillary response 
files that account for this time- varying effect. In general, 
we grouped the spectra into bins of a minimum of 30 
counts, so statistics are appropriat e in the mode l fit- 
ting, which we performed in XSPEC ijArnau dl lT99l . In 
the best-fitting models discussed below, only data from 
0.4 to 8.0 keV are considered, and the inclusion of addi- 
tional model components and their free parameters are 
significant at a minimum of the 95% confidence limit, 
based on an F test. All quoted errors are 90% confi- 
dence for one parameter of interest. 

Observational results from "pure" AGN and "pure" 
starbursts physically motivate the spectral models we 
consider, and we recognize that all components may be 
present here, even in spatially-restricted regions. The 
AGN produces an X-ray continuum that is characterized 
by a power law of photon index T 1.9, a s observed 
in Seyfert 1 galaxies ijNandra fc Pound'3ll994|) . assuming 
the central engines of Seyfert Is and Seyfert 2s are in- 
trinsically identical. With modest column density along 
the line of sight, the soft X-rays will be photoelectri- 

® http:/ /www. astro. psu.edu/uscrs/chartas/xcontdir/xcont. html 
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Fig. 5.— Chandra hard X-ray (4-8 keV) image of NGC 5135, 
smoothed by a Gaussian of FWHM=1'.'5. The image is scaled 
logarithmically. The overlaid intensity contours begin Scr above the 
background level and increase by factors of two. The AGN is the 
strongest source, although resolved emission due to the southern 
starburst concentration is also evident. 



cally absorbed. When the obscuration becomes Conip- 
ton thick, which is the case here, the emergent continuum 
is entirely reflected and is strongly reprocessed, to F « 
(e.g., iKrolik. Madau. & Zvckll994() within the observed 
energy regime. Fluorescent Fe emission is common in 
the X-ray spectra of Seyfert galaxies. The Fe Ka line 
appears at central energy Ec = 6.4 keV in material less 
ionized than Fe XVII. If the iron is nearly fully ionized, 
Ec = 6.7 or 6.9 keV, for He-like or H-like ions, respec- 
tively. Each of these fluorescence lines may be spectrally 
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characterized by a single Gaussian. The starburst pro- 
duces X-rays in both point sources and genuinely diffuse 
gas. Observed starburst spectra tend to be soft, indicat- 
ing multiple thermal components in the diffuse material 
as we ll as a hard continuum (e.g., 'Dah lem et aUlTOQ^ 
l200n() . The X-ray binaries produce the power-law con- 
tinuimi, although the integrated spectrum of the unre- 
solved sources may not be characteristic of any individual 
sources. 

4.1. The AG N Region 

We centered the extraction of the northern central 
source on the Fe Ka center measured in the images. 
We measured the spectrum in a region of only 1'.'2 ra- 
dius to avoid contamination by the southern source. The 
best-fitting model requires two thermal components and 
a Gaussian at soft energies in addition to a flat power 
law continuum and a Gaussian line, which are character- 
istic of the AGN at higher energies. We apply the same 
absorption by a foreground screen to all of these compo- 
nents. Figure El contains the data, model, and ratio of 
data/model, and TableHlists the model parameters. We 
had previously fit the 4-8 keV spectrum l)Levenson et al.l 
12002') with a power law and a single Gaussian to measure 
the Fe Ka line and its underlying continuum. In order to 
retain significant information at energies above th e line, 
we fit the unbinned spectrum using the C-statistic ijGashI 
Il979j) . The Fe line is spectrally unresolved, with equiva- 
lent width EW = 2.4 keV. The power law is flat (F = 0), 
characteristic of a pure reflection spectrum. The large 
EW and flat spectral slope demonstrate that the AGN is 
not viewed at all directly within the Chandra bandpass; 
it must be obscured by Nh > lO^'^cm"^. 
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Fig. 6. — Chandra ACIS-S3 spectrum of the northern central 
source. This region includes the NGC 5135 AGN, whose direct X- 
ray emission is completely obscured. A complex model consisting 
of two thermal components, a power law, and two emission lines fits 
the data well. In the upper panel, the data are plotted as crosses, 
and the model, convolved with the detector response, is plotted as 
a histogram. The ratio of data/model is plotted in the lower panel. 

These results most accurately describe the high- 
energy emission, so we fix these parameters while 
fitting the entire 0.4-8 keV binned spectrum. In the 
complete spectrum, two additional components are 



necessary to account for most of the soft X-ray emis- 
sion. We use th e thermal equilibrium model of Mewe 
and Kaastra j|Mewe^_GTonenschil d, Sz van den OordI 

19851 lArnaud & Rothcnflu'gl IToM 

Mewe. Lemen. fc van den Oord 1986), with updated Fe 
L calculations ((Liedahl 1992) and solar abundance. The 
two thermal components are at temperature kT = 0.05 
and 0.6 keV. We note that the normalization of the 
softer component is several orders of magnitude greater 
than that of the harder component. The normalization 
is proportional to the volume integral of the square 
of the density, and a greatly increased density in the 
low-temperature region is the likely origin of this 
difference. 

Alternatively, we fit the entire 0.4-8 keV spectrum si- 
multaneously. Again, we use the unbinned data and the 
C-statistic to retain some information above the high- 
energy emission line. These parameters are not signifi- 
cantly different from the previous fit and are listed in the 
second row of Tabled The slightly steeper photon index 
(F = 0.3) contributes to the observed 1-2 keV emission, 
and the high-energy line EW is still large (EW — 3.0 
keV). 

Despite the uncertainties in measuring F and its nor- 
malization independently, the Fe EW remains large in 
both the high-energy and full-spectrum fits. Considering 
the errors on these parameters jointly, a steeper photon 
index does require a higher normalization (measured at 
1 keV) , although the magnitude of the high-energy con- 
tinuum change is small. 

The fiat power law of the high-energy model is a 
consequence of reflection. We also considered a more 
complete model based on Monte Carlo simulations of 
Compton reflection by a col d population of electrons 
fMag dziarz fc Zdziarskil 119951 the PEXRAV model in 
XSPEC), which further demonstrates that we observe 
only the reflected continuum and does not improve the 
fits at all. Specifically, we adopted standard values of 
the parameters to which we have no sensitivity, fixing 
the high-energy cutoff Ec — 100 keV and inclination an- 
gle i = 63°. We first considered a fixed intrinsic photon 
index F = 2.0. Fitting the full spectrum we find the un- 
physical result that the magnitude of the reflected spec- 
trum relative to the observed intrinsic spectrum exceeds 
30, and the Fe Ka EW = 2.9 keV. The photon index can- 
not be reasonably constrained. If it is free parameter, we 
find F = 2 (-^3,-2). 

With the spectral resolution and high-energy sensitiv- 
ity of the present data, there is no discernable difference 
between a pure refiection model and the simplified fiat 
power law that we have adopted to describe it. Above 4 
keV, the only formal difference between these two models 
is the Fe absorption edge near 7 keV, which we cannot 
detect. The softer absorption edges of other elements 
in the reflection model are not strong enough to be de- 
tected where thermal emission dominates, despite Chan- 
dra's improved sensitivity at lower energies. 

We considered the possibility that photoionization by 
the AGN might produce the soft X-ray emission in the 
northern region. We attempted to fit the spectrum 
with photoionization models produced with XSTAR^ in- 
stead of the thermal plasma components. We also com- 

^ http:/ /hcasarc. gsfc.nasa.gov/docs/software/xstar/xstar. html 
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Table 1. Spectral Model Paramters 
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Note. — Errors arc 90% confidence limits for one interesting parameter. Parameters tliat are constrained by liard limits are marlicd witli a colon. Fixed parameters are 
marlicd witii f. 

''Column density in units of 10^" cm~^. 
''Temperature of thermal plasma in keV. 

'^Normalization of thermal component in units of 10~^K, where K = (10~^^/(47r_D^)) J nenndV, D is the distance to the source (cm), 
Tie is the electron density (cm~^), and Ufj is the hydrogen density (cm~^). 
"^Normalization of power law in units of 10~* photons keV~^ cm~^ s~^ at 1 keV. 
""Energy of line center in keV. 
'Equivalent width of line in keV. 

^0.5— 2.0 keV model flux in units of 10~^^ erg cm~^ s~^ . 
''2.0—10.0 keV model flux in units of 10"'* ergcm"^ s~'. 

'Parameters and errors on F, A3, and high-energy components of Eii„g and EWune are based on high-energy fit alone. 
■■Parameters and errors are based on full spectral fit, using the C-statistic. 



pared NGC 5135 with the observed ph otoionization spec- 
trum of NGC 1068 l|Kinkhabwala e t al. 2002| because 
the models fail to reproduce some of the line emission 
that is detected in this Seyfert galaxy. None of these pho- 
toionization spectra accurately reproduce the observed 
emission of NGC 5135. In general, NGC 5135 exhibits 
several broad emission peaks that are due to complexes 
of many unresolved lines, characteristic of thermal emis- 
sion. Even at the relatively low spectral resolution of 
these observations, the emission features of all the pho- 
toionization spectra appear too sharp. 

Finally, in addition to the Fe K line at 6.4 keV and the 
thermal emission, the northern source also requires a sec- 
ond line at 1.8 keV. This line is consistent with neutral 
silicon fluorescence, both in central ener gy and equivalent 
width l|Matt. Fabian, fc R.evnolds flQQTl) for the geometry 
that the Fe emission constrains f H5.2l) . We searched for 
other fluorescent lines that are expected to be relatively 
strong (EW > 100 eV). The emission measured near 2.4 
keV is suggestive of neutral sulfur and can be fit with 
an unresolved Gaussian of EW = 340 eV. We caution, 
however, that this component is not statistically signif- 
icant in the model fit, and the central energy is signifi- 
cantly (l.Sfj) higher than the 2.31 keV expected for sul- 
fur. Some of the lower-energy fiuorcscent lines, including 
oxygen, neon, and magnesium, should have large EWs 
with respect to the reflected continuum, but we cannot 
discern them against the strong soft thermal emission. 

The prominent inner-shell fluorescence lines are char- 
acteristic of pure reflection from an optically thick 
medium. Photoionization of optically thin, low den- 
sity material does produce lin e emission, as in the 
X-ray spectrum of NGC 1068 ijBrinkman et alJ I200I 
Hinkhabwala ct al. I2002() . In NGC 5135, however, an 
additional photoionized component would also produce 
many other strong emission lines, which we do not ob- 
serve. As we note above, photoionization alone, in- 
stead of thermal emission, fails to reproduce the soft 
X-ray spectrum. Alternatively, we consider pure coUi- 
sional excitation, which does produce some Si line emis- 
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Fig. 7. — Chandra ACIS-S3 spectrum of the southern central 
source, the strong circumnuclear starburst of NGC 5135. Typical 
of starbursts, the data require a thermal component and a hard 
X-ray continuum. 



sion at s imilar energies , as observed in other starburst 
galaxies llWeave r et al.l l2002') and Compton thick AGNs 
lIBoller et al.ll2003D . The line component we model here, 
however, is required in addition to the contributions from 
the thermal plasmas that are also present and dominate 
the soft X-ray emission. Thus, we conclude that fluo- 
rescent Si is the likely origin of this additional emission 
line. 

4.2. The Compact Starburst 

We extracted the spectrum from a region of 1'.'2 radius 
to isolate the southern central source, which we asso- 
ciate with the strong circumnuclear starburst. The X-ray 
emission is spatially resolved in all the broad bands mea- 
sured above. The best-fitting model requires one thermal 
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component and a power law, both of which are absorbed 
by the same column density. The temperature of the 
thermal component is 0.7 keV, similar to pure star burst 
galaxies. The thermal emission accounts for most of the 
soft X-rays, and the power law contributes all of the X- 
ray emission above 4 keV. The extinction along the line 
of sight, Nh ~ 1-1 X lO'^^ cm~^, is greater than that due 
to the Milky Way Galaxy alone, although it is compara- 
ble to the value Gonzalez Dclgado ct al. ( 1998) measure 
{Nh = 9.9 X 10^" cm^^) at ultraviolet wavelengths in a 
similarly-sized aperture. 

The non-thermal power law is likely the net effect 
of individual sources, primarily X-ray binaries, in the 
starburst. The slope of the power law is very steep: 
r = 2.6. Although this is atypical of any individual X- 
ray binary system, it is characteristic of the integrated 
spectrum of distinct sources measured in face-on star- 
forming galaxies, w hich contain many observable sources 
ijZezas et al.ll2003) . We considered fixing the photon in- 
dex at a value more typical of individual non-thermal 
X-ray sources (F = 1.9), but found the model fit to the 
data significantly worse, even allowing for additional soft 
X-ray components. Figure \7\ shows this compact star- 
burst spectrum with the best-fitting model and ratio of 
data/model, and Tabled lists the best-fitting model pa- 
rameters. 

4.2.1. Abundance Variations 

We also considered thermal plasmas with variable 
abundances. A starburst wind may be enriched if it 
carries the metals that supernovae produce. On the 
other hand, the wind may have low abundances, ei- 
ther because it consists of unenriched halo material or 
because refractory elements are not present in the gas 
phase. Low abundances have been reported in X-ray ob- 
servations of starburst galaxies, but these results may 
be due to the systematic difficulties of measuring multi- 
temperature plasm as or neglecting underlying contin- 
uum contributions llWeaver. Heckman. fc DahlemI EOOJ: 
IStrickland et"ani20?) ir 

We find slight evidence for non-solar abundances in 
the southern region of NGC 5135. Considering variable 
abundances in the thermal component of the best-fitting 
model above (which includes a non-thermal continuum 
contribution), the iron abundance is best fit at [Fe/H] = 
0.39 (+ 0.14, -0.13) [Fc/Hlra, where [Fe/H]© = 4.68x10-5 
(|Anders fc Grevessd Il989t) . We obtain similar results 
if other refractory elements (Na, Al, Ca, and Ni) are 
constrained to have the same relative abundance as Fe: 
[Fe/H] = 0.37 (-f0.17, -0.12) [Fe/H]o. The plasma tem- 
perature remains the same, although the power law be- 
comes somewhat flatter (F = 2.0). Starburst spectra of- 
ten show enhan ced abundances of alpha elements with 
respect to Fe (Martin. Kobulnickv. fc HeckmanI 120021: 



IStrickland et al.l 2003^. Although we cannot simultane- 
ously constrain Fe and alpha abundance variations inde- 
pendently in these data, we do observe this effect gener- 
ally. Above, varying only Fe abundance, Za/Zpe = 2.6 
(-1-1.3, -0.7). Alternatively, fixing Fe at solar abundance 
and allowing the common abundance of O, Ne, Mg, 
Si, and Ca to vary, Za/Zpe = 2.3 (+1.1, -1.0), with 
F = 2.3. The abundance of oxygen most strongly deter- 
mines the value of this ratio, and here we have a dopted 
[O/H]0 = 8.51 X 10-* d Anders fc GrevessellT989l) . 



We also considered multi-temperature thermal plasmas 
with variable abundances without any non-thermal com- 
ponent. These fits were statistically worse and physically 
unlikely. Even with the abundance variations, these fits 
require an extremely high temperature plasma {kT w 10 
keV) to account for the emission at energies above 4 keV. 

4.3. Diffuse Emission 

We extracted two regions of diffuse emission. The first 
(Dl) is the brighter central region of maximum radius 
6'.'5, excluding an ellipse that encompasses the extremely 
bright north and south sources and a point source within 
the region of interest. The second (D2) is elliptical, ex- 
tending 14" X 10" at position angle PA = 118°, excluding 
the central circle of Dl and interior point sources. 
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Fig. 8. — Chandra ACIS-S3 spectrum of central diffuse emission 
(region Dl). A thermal component and a power law are required 
to fit the data. 

We fit Dl with an absorbed thermal component and a 
power law (Figure|Sl Table^l. For the thermal emission, 
kT — 0.63 keV, and the power-law photon index F — 2.2. 
With the total column density the same for all compo- 
nents and a free parameter, we find Nh = 8.4 x 10^° cm~^ 
in the best-fitting model, but it is not significantly differ- 
ent from the Galactic contribution alone, of Nh,mw = 
4.7 X 10^*^ cm~^. The thermal plasma accounts for most 
of the soft X-rays, and the power law is responsible for 
nearly all the hard X-rays. Spectrally, the extended cen- 
tral region is very similar to the southern central source, 
although it has a much lower surface brightness. The 
ratio of hard to soft luminosity is somewhat larger in 
Dl, compared with the southern source. In this case, 
normal disk sources, in addition to those in the areas 
of significant star formation, likely become increasingly 
important. 

We considered abundance variations in the Dl spec- 
trum but found no evidence for non-solar abundances. 
Similar to the results of the South region, excluding the 
non-thermal component cannot reproduce the higher- 
energy spectrum of Dl, even using multi-temperature, 
variable abundance models. 

The spectrum of the largest region, D2, also requires 
a single-temperature plasma. This component is much 
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Fig. 9. — Chandra ACIS-S3 spectrum of very extended diffuse 
emission (region D2). A single thermal component and power law 
fit the data. 



cooler {kT = 0.25 keV) and more absorbed {Nh = 4.2 x 
10^^cm~^). The required second component is not well 
constrained. We model it as a power-law, fixing F — 2.2, 
the best-fitting value of F in the Dl spectrum. The best- 
fitting model of D2 is plotted with the spectrum in Figure 
1^ and the model parameters are listed in Tabled 

A multi-temperature medium may be a more accurate 
description of physical conditions of this region, possibly 
accounting for the unphysical conclusion that this region 
is more obscured than Dl. With the limited sensitivity 
of these data, however, more complex models are not sta- 
tistically significant. Examining a model that includes a 
second thermal plasma component, we find kT2 ~ 0.9 
keV. With this higher-temperature component account- 
ing for the iJ « 1 keV emission, the normalization of 
the very soft component is lower. The column density 
is also slightly lower {Nh = 3.1 x lO^^cm"^), although 
still larger than that measured in Dl. 

4.4. Additional Sources 

We detect several other weaker sources within 2' of 
the center of NGC 5135, which we list in Table |21 in 
order of distance from the nucleus. The last three are 
unlikely to reside within the galaxy, but we include them 
for completeness, since they are located within a few arc- 
seconds of the optically-detected spiral arms. These ad- 
ditional sources are interesting as possible intermediate- 
luminosity X-ray objects (IXOs), those with X-ray lumi- 
nosity Lx > lO'^^ergs"^. Exceeding the Eddington lu- 
minosity of a stellar-mass black hole, they could be single 
sources if they are more massive, accrete above the Ed- 
dington rate, or emit anisotropically. Although IXOs are 
obser ved in gal axies without intense ongoing star forma- 
tion ijColbert fc Ptak .2002) . in galaxies having high star 
formation rates, a significant popul ation of IXOs is corre - 
lated with the star formation rate iColbert et al.ll200^ . 
Not all of the sources listed in Table |21 are intermediate- 
mass {10^~10^Mq) black hole candidates, however. Four 
of the NGC 5135 sources (2, 3, 5, and 6) are certainly 
spatially extended, and even the unresolved sources en- 



compass spatial scales of several hundred parsecs, so they 
could be clusters of multiple sources. 

The net source counts in the total (0.3-8 keV), soft 
(0.3-2 keV), and hard (2-8 keV) bands are fisted in 
columns 4 through 6, respectively, where the local back- 
ground measured in a surrounding annulus has been sub- 
tracted. We use these energy bands, rather than the 0.5- 
2 and 2-10 keV that are more common in the literature 
and presented elsewhere in this work to avoid errors in 
the model-dependent conversion of the observed data to 
unmeasured energies. In some cases, the total number of 
background-subtracted counts detected does not equal 
the sum of the soft- and hard-band measurements be- 
cause of rounding errors. None of the sources is bright 
enough to distinguish the true physical character of its 
emission, but soft and hard count rates indicate varia- 
tions of spectral hardness in the different sources. We 
note that many of these sources are extremely soft, un- 
like the harder emission typically measured in IXOs. The 
emission of these soft sources could therefore be thermal, 
due to concentrations of supernovae and stellar outflows, 
although we cannot rule o ut accretion, which produc es 
soft spectra in some IXOs lITerashima fc Wilsonll2?)03l) . 

To estimate fluxes and luminosities consistently, we 
model each spectrum as an absorbed power law, even 
though this may not be physically appropriate. We fix 
the photon index F = 1.8, and allow absorption to vary 
in excess of the Galactic column density along this line 
of sight. The observed fluxes and luminosities (not cor- 
rected for absorption) in the 0.3-8 keV band are listed 
in columns 7 and 8 of Table 13 where the luminosity is 
calculated for the 59 Mpc distance to NGC 5135. 

5. PHYSICAL INTERPRETATION 

5.1. Identification of X-ray Emission Sources 

The X-ray spectral and spatial information together 
identify the northern source as the AGN of NGC 5135. 
The most compelling evidence is the Fe Ka line, which 
stellar processes cannot produce with such a large lu- 
minosity or equivalent width. The spectral aperture by 
necessity covers a relatively large area, so in addition to 
the AGN, it includes other emission sources. Therefore, 
the energy-restricted Fe-line image is required to spa- 
tially constrain the AGN. In this uncontaminated image 
of the central engine, we demonstrate that the AGN is 
unresolved. 

At higher spatial resolution, we also identify the AGN 
located north of extended stellar emission in Hubble 
Space Telescope (HST) optical and near- infrared (NIR) 
images. NGC 5135 was observed for 500s with the Wide- 
Field and Planetary Camera 2 (WFPC2) through the 
F606W filter on 1995 February 15 and for 320s with 
NICMOS Camera 2 through the F160W filter on 1998 
May 28. In both the optical and NIR, the nucleus is 
located at the apparent center of relatively bright emis- 
sion that extends over scales of 5". In the latter, at 
1.6/im, it is unresolved and the brightest source. In the 
former, the nuclear emission is resolved, including a sig- 
nificant contribution that extends over physical scales 
greater than 10 pc. The strongly-obscured AGN is not 
the brightest source observed within the 606W filter 
bandpass. In the ultraviolet, the AGN is a weak source 
compared with nearby stellar emission, as observed at 
2150A with the F210M filter of HSTs Faint Object 
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Table 2. Additional Sources 



Source 


R.A. 


Decl. 


Counts 


Counts 


Counts 




Lo.3- 




(J2000) 


{J2000) 


(0.3-8 keV) 


(0.3-2 keV) 


(2-8 kcV) 


(lO^-"-^ crgcm"^ s^-"-) 


(10^^ erg 


1 


13 25 43.7 


-29 49 56 


46 


43 


3 


7.7 


3.2 


2 


13 25 44.7 


-29 50 04 


41 


26 


15 


7.5 


3.2 


3 


13 25 44.8 


-29 50 14 


25 


20 


4 


6.0 


2.5 


4 


13 25 42.6 


-29 50 06 


13 


13 





6.4 


2.7 


5 


13 25 45.3 


-29 50 08 


20 


14 


5 


4.9 


2.1 


6 


13 25 45.2 


-29 50 12 


8 


7 





3.6 


1.5 


7 


13 25 44.1 


-29 49 33 


8 


7 


1 


2.7 


1.1 


8 


13 25 45.3 


-29 49 35 


6 


5 


1 


1.9 


0.80 


9 


13 25 48.0 


-29 49 48 


35 


31 


5 


11 


4.6 


10 


13 25 46.7 


-29 51 43 


9 


7 


2 


0.98 


0.41 


11 


13 25 44.1 


-29 48 10 


68 


46 


22 


13 


5.5 


12 


13 25 52.8 


-29 50 09 


30 


25 


10 


5.4 


2.3 



Note. — Units of right ascension are hours, minutes, and seconds, and units of decHnation are degrees, arcminutes, and arcseconds. 



Camera for 1963s on 1995 July 17. (We note that the 
AGN is located approximately 0'.'2 west and 1'.'7 north 
of the UV-brightest point and is not the location that 
iGonzalez Deleado et al.l l)1998j) tentatively identified as 
the nucleus. In fact, their UV spectral aperture proba- 
bly excluded most of the AGN light.) 

Because of errors in the absolute astrometry of both 
Chandra and HST, the location of the northern X-ray 
source does not appear to be exactly coincident with the 
UV, optical, or NIR AGN, based on the telescopes' re- 
ported astrometry alone. In Figure [TUI we illustrate the 
approximate relationship between the X-ray and optical 
emission, having shifted the X-ray contours by 1'.'8 to 
align the northern source with the optical AGN and re- 
sampled onto a common pixel scale. A further clockwise 
rotation of the X-ray image may also be suggested, but 
the simple planar shift alone illustrates the strong stellar 
contribution south of the AGN, near the southern X-ray 
source. 

The AGN is not the only source of X-ray emission 
within the 340-pc northern spectral aperture. Mass is 
strongly concentrated in the nuclear region, so it is a lo- 
cation of the most intense star formation. At the higher 
resolution of WFPC2, we distinguish many individual 
stellar sources at optical wavelengths on sub-arcsecond 
scales. Stellar processes, including supernovae, stellar 
winds, and large-scale outflows, are likely responsible for 
the soft thermal X-ray emission, even in the immediate 
surroundings of the AGN. Most of the thermal emission 
in the northern region is due to the 0.6-keV component, 
which is similar to the only significant thermal compo- 
nents of the South and Dl regions. Thus, we suggest that 
even within the northern aperture, the thermal emission 
is stellar, not due to the AGN. 

At 6 cm, VLA observations show concentrated emis- 
sion at a = 13^25™43."94, 5 = -29°50'02'.'8 (J2000), 
with s ome low-level emission ext ending toward the north- 
east ljUlvestad fc WilsonI [T989j) . The radio core ap- 
pears to be offset 2'.'9 south of the AGN, which is 
located in the continuum-subtracted Fe Ka image at 
a = 13''25'°43."97, 5 = -29°50'00'.'0 (J2000). The ra- 
dio core appears aligned instead with the southern X-ray 
source (within 0'.'3), suggesting that the concentrated ra- 
dio emission is due to the starburst. 



T3 




4 2 -2 -4 

arcseconds 



Fig. 10.— HST optical (6O6OA) image with C/iartdra total-band 
contours overlaid, scaled linearly. North is up and east is to the left. 
The AGN is located at the origin, with the X-ray emission shifted 
to align at the northern X-ray source. In both images, the AGN is 
north of prominent, extended emission. The high-resolution optical 
data also illustrate that stellar sources, in addition to the AGN, 
are contained within the small nuclear X-ray aperture. 



The southern source is simply a starburst. The ther- 
mal emission and the hard continuum are typical spectral 
components measured in local starbursts. The power-law 
continuum in the starburst region is therefore unrelated 
to the nearby AGN in this galaxy. Although only one 
thermal component is significant in this region, the real- 
ity is likely more complicated, with a range of tempera- 
tures actually present (Strickland & Stevens 2000). 

The starburst is viewed approximately face-on. The 
southern X-ray extraction includes most of the high 
surface brightness region measured in the UV con- 
tinuum ijGonzalez Delgado et al1ll998D . Edg e-on star- 
bursts, such as M82 and NGC 253, indicate the 
three-dimensional stru cture that the aperture covers 
ijStrickland et al.l2003l ). Specifically, the X-ray- and UV- 
brightest regions of intense star formation are located 



Accretion and Outflow in NGC 5135 



9 



at the base of an outflowing superwind, which extends 
conically out of the plane of the galaxy. The 0.7-keV 
thermal component is most similar to the diffuse outflow 
measured in the edge-on cases, where it can be spatially 
isolated. The 200-pc radial scale of NGC 5135's bright 
starburst is comparable to the intense cores of edge-on 
starbursts, so we expect that the outflow extends on sim- 
ilar physical scales — on the order of a kpc or two — along 
the line of sight. 

The diffuse region Dl is spectrally similar to the in- 
tense starburst, with the same hard power-law and soft 
thermal components, although at greatly diminished sur- 
face brightness. These two regions may be genuinely re- 
lated, representing the outer walls of the expanding con- 
ical outflow, which is projected onto the starburst base 
at small scale height and has a larger projected radius 
at large scale height. The more extended component, 
D2, arises in the larger (4 kpc scale) star-forming area 
that the optical and UV HST images show, but X-ray 
emission does not cover the entire galaxy disk. Thus, 
even this extended component is closely tied to the sites 
of significant star formation, so this diffuse emission is 
most likely a halo rather than disk component of NGC 
5135's interstellar medium. 

5.2. Fluorescent Iron Emission 

The properties of the Fe line measured in the AGN 
spectrum are consistent with fluorescence in a cold, "neu- 
tral" medium. The spectrally unresolved line limits the 
temperature to T < 2 x 10^ K. The energy of Fe Ka is 
about 6.4 keV for all ions up to Fe XVII, and only reaches 
6.7 and 6.9 keV for He-like and H-like Fe, respectively. 
Thus, the Fe line we observe in NGC 5135 is not broad or 
ionized enough to arise in a scattering medium that views 
the intrinsic continuum directly. In these more sensitive 
and higher resolution data, we do not find any evidence 
for th e strong ionized Fe lines reported bv lTurner et al.l 
l|1997t) based on ASCA observations. The hard X-ray 
fluxes of these two observations separated by over 5 years 
are similar, so AGN variability is an unlikely explanation 
for this discrepancy. 

The large fluorescent iron equivalent width, EW = 2.4 
keV, further constrains the geometry of material close 
to the AGN. In the simplest case, where both the con- 
tinuum and Fe lines are viewed directly, EW < 250 eV 
((Krolik & Kallman 19 87), which is ty pical of Seyfert 1 
galaxies (e.g., iNandra fc Poundsll994j) . Reaching EW > 
1 keV requires high col umn density {Nh > 10^"^cm~^) in 
the fluoresci ng reg ion ijGhisellini. Haardt. fc Mat'3ll994l: 
IKrolik et alJ[l99l . 

Detection of an extremely large EW strongly con- 
strains the geometry of th e obscu rer. Extending t he sim- 
ulations of IKrolik etHI l)1994|) . iLevenson etldl l|2002fl 
show that variations of geometry, not iron abundance 
or intrinsic spectral shape, are required to produce EWs 
signiflcantly greater than 1 keV. Specifically, the largest 
EWs require that the solid angle covering fraction of the 
AGN is also large. Equivalently, we may describe this 
situation as a small opening angle of a toroidal obscur- 
ing region. In NGC 5135, for example, the measured EW 
constrains the opening angle, 9 < 20°, or covering frac- 
tion > 90%. Figure ITTl illustrates the effect of variations 
of 9 and Thomson optical depth along the line of sight, 
r, on the EW. The Fe Ka EW measured in NGC 5135 
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Fig. 11. — Fe Ka line EW as a function of Thomson opti- 
cal d epth and opening ang le of the obscuring region ( Krolik ct all 
[T99l'Leven son et aH200a) . The extremely large EW of NGC 5135 
(dashed line) requires an opening angle 9 < 20°, or equivalently a 
covering fraction > 90%. 



is also plotted. 

The large covering fraction may be related to the cir- 
cumnuclear starburst in this case. The mechanical en- 
ergy of stellar winds and supernovae may help to "puff" 
the obscuring torus of material to become more spherical 
than disk-like. Whatever the origin, the large covering 
fraction means that very few lines of sight to this galaxy 
have a direct view of the central engine. Galaxies simi- 
lar to NGC 5135 are therefore extremely unlikely to be 
observed as Seyfert Is. 

The luminosity of the Ka line is related to the intrinsic 
luminosity of the AGN, although the uncertainty in the 
exact correlation is large. Based on the Ka luminosity, 
we estimate the intrinsic 2-10 keV luminosity of the AGN 
L2-io,int = 1 X lO^'^ergs"^, considering o pening angles 
9 < 20° in the Monte Carl o simu lations of IKrolik et alJ 
(1994) and ILevenson et "all l)2002|) . As expected in this 
Compton thick AGN, less than 1% of the intrinsic X- 
ray luminosity emerges from the obscured region at 
hard X-ray energies; the observed 2-10 keV luminosity, 
L2-10, AGN = 3 X 10''°ergs~^. Scaling the predicted in- 
trinsic hard X-ray luminosity to the obser ved spectral 
energy distributions of unobscured quasars (jElvis et al. I 
1994), we estimate the bolometric luminosity of the AGN 
below 10 keV; LboLAGN = 3 x 10'''' ergs~^, which is only 
half of the bolometric luminosity of the galaxy as a whole 

5.3. Absorption and Geometry 

Along the direct line of sight to the AGN we find 
Nh > lO^^cm"^, based on the large Fe Ka EW and 
the absence of a continuum characteristic of a directly- 
viewed AGN, even at hard X-ray energies. The optical 
[O III] A5007 luminosity relative to the 2-10 keV emission 
offers further support for the Compton thick classiflca- 
tion of NGC 5135. The [O III] emission is expected to 
arise outside the obscuring region, while the emergent X- 
ray continuum may be suppressed along the line of sight. 
A ratio L2~io/ L\o jjj] < 1 is chara cteristic of Comp- 
ton thick cases l|Maiolino et aljri998(l . In a 2-4" aper- 
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Table 3. Region Luminosities 



Region 


io.5-2 


^2-10 




(10''0 ergs-1) 


(10*"crgs-i) 


NGC 5135 (total) 


18 


19 


North (AGN) 


2.6 


9.4 


South (SB) 


4.6 


1.5 


Dl 


3.1 


1.6 


D2 


1.5 


0.5 



ture, iWhittld l)1992|) measures a flux /[o ui] = 2.2 x 
10"^'^ ergcm"^ s~^. Cor recting: for reddening based o n 
the Balmer decrement ijGonzalez Delgado et alJll998D . 
we find L2-io/^[o ///] = 0.08. Starburst emission cer- 
tainly contributes to the measured [O III] luminosity, but 
this result is at least consistent with Compton thick ob- 
scuration of the AGN. 

The starburst itself is responsible for some of this ob- 
scuration. In general, the circumnuclear starburst re- 
quires a concentration of mass, some of which is lo- 
cated in front of the active nucleus. In this case, we use 
the star formation rate measured within the 1'.'7 x 1'.'7 
aperture of the Goddard High Resolution Spectrograph 
(GHRS) on HST to estimate the column density near 
the nucleus. A star- formation rate SFR w 2Moyr~^ 
({Gonzalez Delgado et al.llT998() within this aperture cor- 
responds to a mean gas c olumn dens i ty Nu ~ 10^^ cm~^, 
using the conversion of iKennicuttI |1998). In extreme 
star-forming galaxies, such as Arp 220, the observed 
SFR alone produces Compton thick column densities 
along the line of sight. 

We find slightly lower column densities in the Chan- 
dra observations of the starburst areas, including the soft 
thermal emission in the North region, because we do not 
directly probe the line of sight to the very center of NGC 
5135. We do, however, always measure absorption in ex- 
cess of the Galactic column density alone. The X-ray 
measured obscuration, Nh = 4 x 10^^ cm^^, is greater 
than that measured within the starburst itself in the UV, 
where E{B ~ V) ^ 0.20, or Nr = 1.2 x lO^^ cm^^ in 
the comparable GHRS aperture. The UV extinction re- 
lies on modelling the observed emission features assum- 
ing a foreground screen, and therefore represents only 
a lower limit on the total extinction through the en- 
tire starburst region. Sources located at r < 1 domi- 
nate the observed emission. Thus, it is not surprising 
that the detected AGN, located behind the entire ob- 
scuring starburst, is even more absorbed than the ther- 
mal emission alone, or that the UV measurement yields a 
smaller column density. The UV-determined extinction 
is comparable to the obscuration measured in the X-ray 
spectrum of region Dl. The column density of the D2 
extraction is greater. Wc note that in other UV spec- 
tra, [Gpnzalcz Delgado ct al. (1998) also measure greater 
obscuration in larger apertures in NGC 5135 and other 
Seyfert 2/starburst galaxies. 

6. FAR-INFRARED RELATIONS AND SPECTRAL ENERGY 
DISTRIBUTIONS 

NGC 5135 is classified as a normal Seyfert 2 galaxy 
on the basis of optical emission-line ratios, yet at most 
wavelengths, the AGN contributes only a small fraction 



of the detectable luminosity. Even in X-rays, the star- 
burst is significant. Table 01 lists the soft and hard X-ray 
luminosities of the four spectroscopic regions and a 1' 
aperture, which represents the total emission of NGC 
5135. 

The far-infrared (FIR) emission of NGC 5135 is char- 
acteristic of star-forming galaxies, both in spectral shape 
and relative magnitude. While active galaxies typi- 
cally exhibit 25- to 60-Mm flu x ratios /25//6O > 0.26 
(jde Griip. Lub. fc MilevI 1198 7) in the arcminute-scale 
apertures of the Infrared Astronomical Satellite {IRAS), 
F25/F60 = 0.14 in NGC 5135. The FIR luminosity 
LpiR — 5.5 X lO^'' ergs~^, and the total 8-1000/im lumi- 
nosity L tr = 7.1 X lO'^'^ ergs~^, following the pre scrip- 
tions o f 'Helou. Soifer. fc R owan- RobinsonI (|1985') and 
iSanders fc Mirabel (,1996) to calculate these luminosities 
from IRAS fluxes. For starburst galaxies, Ltoi ~ Lir, 
and in this case, the observed IR luminosity is more than 
twice the predicted intrinsic luminosity of the buried 
AGN. 

The bolometric luminosity is related to the star 
formation rate, with SFR = L^oil"^-^ x lO''^ 
(|Leitherer fc HeckmanI I1995D . Here we find SFR = 
25MQyr~^. As expected, this is greater than the UV 
SFR measured within the central 500-pc GHRS aper- 
ture. It is less than the extinction-corrected SFR « 
50Mq yr""'^ derived from the International U ltraviolet 
Explorer {lUE) data ({Gonzalez Delgado et al.ll998^1 . but 
because the correction for extinction is highly uncertain 
and large, UV- and IR-mcasurcd SFRs often disagree 
{Mcurcr, Hcckman, fc Calzctti 1999). 

The ratios of total X-ray to FIR luminosity measured 
in NGC 5135 are characteristic of starburst galaxies, not 
AGN. In Seyfert 1 and Seyfert 2 galaxies that do not 
contain starbursts, L2-10/ Lpjn > 10~^, while this ra- 
tio tends to be much lower in starburst and ultralumi- 
nous infrared galaxies ()Levenson et al . 2001a: ?tak et alJ 
IMi. We find L2-w,totai/LFiR = 4.8 x 10-^ in NGC 
5135, where £2-10,(0*0/ is the total observed 2-10 keV lu- 
minosity. Considering only the non-AGN contribution, 
which we take to be the North region contribution sub- 
tracted from the total, L2-^Vl.^^Bl L pjR = 2.4 x 10~^, ex- 
actly the linear relationship IRBualli. Comastri. fc SettH 
(j2003) fit to a local sample of starburst galaxies. The 
soft X-ray/FIR ratios of NGC 5135 iLo.5-2,totai/ Lfir. = 
4.6 X 10"'' and L^,^-2.sb/ Lfir^ — 3.9 x 10~^) are slightly 
higher than the starburst average, but within the typical 
scatter of the data. 

Figure [T^ illustrates the spectral energy distributions 
(SEDs) of extended emission and the AGN separately. 
Across the spectrum, NGC 5135 as a whole (plot- 
ted with open circles) appears most si milar to dusty 
starburst galaxies, whose average SED (jSchmitt et al. I 
Il997r is plotted, scaled to the luminosity of NGC 5135. 
Here we estimate the average X-ray SED from the FIR 
measu rements using the relationship of Ranalli et aQ 
()200 ^. The large-aperture data are compiled from the 
I RAS Faint Source Catalog, ground-based observations 
oflGlass fc MoorwoodI (119851) ■ and /KE meas urements of 
iStorchi-Bergmann. Kinnev. fc Challi^ ()1995j) . The "ex- 
tended" X-ray data that are indicated at low spectral 
resolution come from the measurements of the entire 
galaxy excluding the central 2'.'5. (This excluded region. 
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Fig. 12. — Spectral energy distribution of NGC 5135, showing 
large-scale emission (open circles) and the spatially isolated AGN 
contributions (filled diamonds). Thermal contributions to the nu- 
clear X-ray spectrum are marked separately (open diamonds). In 
general, the galaxy appears very similar to dusty starburst galax- 
ies (dotted line;|Schmitt ct al. 1997). For comparison, the average 
radio-quiet quasar spectrum I E lvis et al. 111994) scaled to the in- 
trinsic bolometric luminosity of NGC 5135 is shown (solid line). 
Simple foreground extinction of this spectrum cannot reproduce 
the unresolved emission. The emergent X-ray spectrum is strongly 
absorbed and reprocessed, and the SED near optical wavelengths 
indicates the presence of some additional foreground extinction. 



slightly larger than the AGN aperture, was selected to 
avoid all AGN contamination.) The solid filled X-ray 
points are measured in the North region and include only 
the spectral components characteristic of the AGN. The 
data points are derived from monochromatic fluxes. The 
highest energy data plotted (at 5 and 8 keV) indicate 
continuum flux only and exclude the Fe Ka line. The 
open diamonds represent the thermal emission measured 
in the nuclear region, which we do not attribute imme- 
diately to the AGN. 

In the UV, optical, and NIR, we measure the AGN 
alone by fitting the appropriate two-dimensional point- 
spread function (PSF) to images from the HST obser- 
vations mentioned above. We modeled the optical and 
NIR PSFs with Tiny Tim^ and used the observation of 
a point source for an empirical UV PSF. The formally 
best-fitting PSF subtraction within a limited area pro- 
duces unphysical residuals that decline toward the cen- 
ter and th erefore yields an upper limi t on the unresolved 
emission (Ridg wav fc Stockt on 19971. We constrain the 
residuals (which represent genuinely extended emission) 
to be monotonic across small interior regions for a more 
realistic estimate of the AGN contribution. This con- 
straint reduces the AGN flux by about 10% in the NIR 
and 25% in the optical. It is most restrictive in the UV, 
where the truly unresolved emission appears to be only 
one-third of the flux within the central region. For the 
unresolved emission, we find 6.7 x 10"'^, 9.0 x 10~^, and 
1.1 X 10^^ Jy in the NIR, optical, and UV, respectively. 

Where the AGN can be spatially isolated, its spec- 
tral energy distribut ion does not matc h that of observed 
radio-quiet quasars ijElvis et al. II199H) . A simple screen 

® http:/ /www. stsci.edu/software/tinytim 



of obscuration in the foreground of the average unob- 
scured quasar also fails to describe the AGN SED. The 
Compton thick reprocessing, which strongly alters the 
X-ray spectrum, is essential. Although spatially unre- 
solved, the observed lower-energy emission must also be 
reflected; we would not detect the AGN directly at these 
wavelengths behind a column density of 10^'*cm~^. The 
steep decline of the relative luminosity from the NIR 
through UV, however, indicates that some additional 
foreground obscuration is also present. 

This example illustrates the difficulty of identifying an 
AGN in the presence of a starburst, especially more dis- 
tant cases where spectral apertures cover larger phys- 
ical scales and signal-to-noise ratios are diminished. 
The characteristic optical emission lines are easily di- 
luted when off-n uclear sources account for a la rger frac- 
tion of the light J Moran. Filippenko. fc Chorn ock 200^ 
iKauffmann et all 20031) . or when a strong~starburst is 
also present l)Levenson et al.ll2001bj) . 

Only above 4 keV does the AGN dominate the de- 
tectable emission from NGC 5135, yet simple modeling of 
the total X-ray spectrum would yield misleading results. 
The detected 2-10 keV luminosity could be produced by 
a starburst alone. The spatially-integrated spectrum is 
relatively flat from 0.4-8 keV, offering no indication of 
the distinct physical origin of the soft and hard compo- 
nents. The strong obscuration of the AGN is not at all 
evident in the spatially-integrated spectrum. A crude 
technique such as spectral hardness ratios would incor- 
rectly suggest the presence of a weak yet unobscured 
AGN. The example of NGC 5135 demonstrates how such 
Compton thick AGN may therefore remain successfully 
hidden at most wavelengths where surveys have been 
performed, particularly when they are located at large 
redshift. 

Compton thick examples are truly normal AGN, 
comprising roughly half of samples that are not 
selected on the basis o f their X-ray emission 
(jRisahti. Maiolino. fc Salvatil Il999j) . They do not 
necessarily appear to be X-ray bright below 10 keV. 
Stellar processes further complicate their soft X-ray 
spectra, since the large reservoirs of gas that block the 
AGN can also form stars, which subsequently produce 
soft thermal emission. 

Though difficult to identify, such obscured sources 
may contribute significantly to the cosmic X-ray back- 
ground. Synthesis models that combine a distribu- 
tion in both absorption to the central engine and red- 
shift iSetti fc Woltieill989l:lMadau. Ghisellini. fc Fabiai] 
119941 IComastri et all 119951) are essential to reproduce 
the background spectrum, which peaks around 30 keV. 
While discrete sources account for 80% or more of the 
XRB at energies up to 10 keV llBrandt et alJ 1200 It 
IGiacconi etlll I200H IHasinger et alJ \2001\i . these ob- 
served redshift and column density distributions do not 
spectra fly match t he X-ray background at higher en- 
ergies (^Gim"200^. Further more, theoretical models 
fPabian et al. ..1998,') suggest that up to 90% of the ac- 
cretion power in the Universe is hidden from direct view. 
NGC 5135 offers a local demonstration of how such 
sources may remain effectively buried behind large col- 
umn densities and along most lines of sight. 

7. CONCLUSIONS 
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In analyzing this Chandra ACIS observation of the 
Seyfert 2/starburst galaxy NGC 5135, we spatially and 
spectrally distinguish both the AGN and the stellar ori- 
gin of its X-ray emission. Images reveal two strong and 
concentrated central sources. The northern of these is 
the AGN, while the southern is associated with star for- 
mation and likely represents the base of an outflowing 
starburst superwind. Recognizing the contributions of 
both the AGN and stellar processes, we apply physically- 
motivated models to spectroscopy of several distinct re- 
gions of interest. In addition to the two central sources, 
these regions include two extended areas, which likely 
represent a hot halo component of NGC 5135. In the nu- 
clear region, the AGN produces a very flat (F = 0) power 
law continuum and extremely prominent (EW = 2.4 keV) 
Fe Ka fluorescence line. These spectral characteristics 
demonstrate that the central engine is highly obscured, 
behind Nh > lO^'^cm"^. The intrinsic emission emerges 
only after being strongly reprocessed and is diminished 
by factors of 100 as observed in the Chandra bandpass. 
This central region encompasses a physical scale of 340 
pc, which includes sites of prominent star formation that 
are viewed directly in UV through NIR wavelengths. 
These stars are responsible for the soft thermal emission 
we measure in this region. 

The remainder of the X-ray emission, including the 
southern source, is spatially resolved. We attribute this 
emission to star formation. Prominent soft thermal X- 
rays characterize these spectra. They also require some 
additional continuum components, which are responsible 
for their higher-energy (E > 4 keV) luminosity. Unre- 
solved populations of X-ray binaries likely produce this 
emission. 

The vast majority of the detected soft X-ray luminosity 
from NGC 5135 (Lo.5-2 = 1-8 x lO'^^eTgs~^) is spatially 
extended. The southern source accounts for fully one- 
quarter of this emission. The AGN itself accounts for 
almost none of it, given that nearly all of the soft X-rays 
within the nuclear aperture are thermal. The AGN is 
directly responsible for roughly half the detected hard 
X-ray luminosity (L2-10 ~ 1-9 x lO^^ergs"^). Based 
on the luminosity of the Fe Ka line, we estimate the 
intrinsic bolometric luminosity of the AGN: Lhoi,AGN = 
3 X 10'*^ erg s^^. This AGN luminosity is only half of the 
bolometric luminosity of the galaxy as a whole, Li,oi ~ 
7 X lO'^^ergs"^, based on the IR luminosity. 

We suggest that the starburst is related to the obscu- 
ration of the AGN. The starburst itself directly hides the 
AGN, at least in part. The measured star formation rate 
alone can account for 10% of the indicated absorption. A 
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starburst requires large reservoirs of gas in the center of 
the galaxy, and the dynamic conditions that concentrate 
this material also serve to further obscure the AGN. The 
starburst may also be related to the activity of the central 
engine, with the instabilities that lead to star formation 
aiding accretion. 

Although optically classifled as an ordinary Seyfert 2 
galaxy, the Fe Ka emission is the only X-ray feature that 
truly identifies the presence of an AGN. The extremely 
large EW of the line requires a very large (90%) covering 
fraction of material, allowing few direct views of the cen- 
tral engine from any line of sight. The Compton thick 
obscuration severely diminishes the hard X-ray emission 
that is a useful discriminant of less-obscured AGN. 

NGC 5135 represents a local case study of some of the 
complex reality that is also likely present in more distant 
examples. Even in X-rays, stellar processes are signifi- 
cant, and combined with Compton thick obscuration, a 
deeply buried AGN will not be obvious in data that have 
low effective spatial or spectral resolution. These are 
particular difficulties in the study of ultraluminous in- 
frared galaxies, in which star formation is known to be 
important, and also in the resolution of the X-ray back- 
ground, where large obscuration along most lines of sight 
diminishes the likelihood of identifying the hidden AGN 
that are essential to produce the background spectrum 
at higher energies. 
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